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Abstract 

The Kondo problem of two impurities in ID strongly correlated electron system 
within the framework of the open boundary Hubbard chain is solved and the impu- 
rities, coupled to the ends of the electron system, are introduced by their scattering 
matrices with electrons so that the boundary matrices satisfy the reflecting integra- 
bility condition. The finite size correction of the ground state energy is obtained 
due to the impurities. Exact expressions for the low temperature specific heat con- 
tributed by the charge and spin parts of the magnetic impurities are derived. The 
Pauli susceptibility and the Kondo temperature are given explicitly. The Kondo 
temperature is inversely proportional to the density of electrons. 
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In recent years the Kondo problem J5J in the low dimensional strongly corrected electron systems 
has been the focus of intense activity and the low temperature properties of such systems in one 
spatial dimension can be described in terms of a Luttinger liquid ||. The electron- electron 
interaction is discarded in the original treatments of Kondo problem, which is reasonable in three 
dimensions where the interacting electron system is described by Fermi liquid. The recent advances 
in semiconductor technology enable to fabricate very narrow quantum wire which can be considered 
as one-dimensional system and furnishes a real Luttinger liquid. Edge states in a 2D electron gas 
for fractional quantum Hall effect can also be regarded as Luttinger liquid ||. Around these 
subjects there are intense efforts and much progress has been made from different approaches by 
the use of the bosonization and renormalization techniques [Q(5), boundary conformal field theory 
|(|], the scaling method J?J and the density matrix renormalization group calculations. 

Kondo impurities play an important role in strongly electron models, especially in one dimen- 
sional system, which is due to only a small amount of defects may change the properties of the 
system, and the impurities usually destroy the integrability of the pure system when they are 
introduced. In the pioneered work, Andrei and Johannesson introduced a spin-5* impurity into 
the isotropic Heisenberg chain with preserving the integrability of the system. It was developed in 
Ref. [9] by Lee and Schlottmann for the host chains with the higher spin values. Eckle, Punnoose 
and Romer studied some physical properties of the model [|l0| . Recently, Bediifig et al has solved 
an integrable impurity model where the impurity is coupled to t — J chain with both spin and 
charge degrees of freedom . The impurity is introduced through a local vertices as in Ref. [8] . 
Schlottmann and Zvyagin introduce the impurity in supersymmetric t — J model via its scattering 
matrix with the itinerant electrons |l2| and the Hamiltonian of the system can be constructed by 
the transfer matrix. They have discussed also the magnetic impurities embedded in the Hubbard 
model jl3| . For all these cases the periodic conditions have been imposed on the electron host and 
the impurity models suffer the lack of backward scattering. Generally speaking, the impurities 
can divide a 'pure' system into several part with the impurities located at the ends of the every 
part. And based on Kane and Fisher's observation ||, we see it is advantageous to use open 
boundary problem with the impurities at open ends to study the problem of impurities coupled 
with strongly-correlated electron system. Indeed, this program has been used for S— interacting 
fermi system in Ref. [14] , t — J model in Refs. [15-18] and Heisenberg spin chains in Ref. |l9|-^l]j 
with the open boundary conditions. 

The one-dimensional Hubbard chain describes the strongly correlated electron system and there 
are many works since the exact solution was found by Lieb and Wu p2| with the periodic condition. 
The one-dimensional Hubbard chain with open boundary conditions jp3|, of boundary chemical 
potentials |^4| and boundary magnetic fields ^,^6| , is solvable also by the use of the Bethe ansatz 
method. Recently, E/31er and Frahm (27j study the singularities in X-ray absorption spectra of 
the one-dimensional Hubbard model and the t — J chain and some thermodynamic properties of 
the model at zero temperature. There are various complex solutions based on the Bethe ansatz 
equations and described the boundary bound states for antiholons, spinnons, and pairs of electrons, 
respectively pq] . The finite-size spectrum for the relevant boundary conditions is given in Refs. 
[29,24,27]. The effect of boundary point contact potentials on the coherent mesoscopic oscillations 
is studied by Frahm and Zvyagin |3^] for the strongly correlated quantum wires. 

In this letter, we devote to study the magnetic impurities in the open boundary Hubbard chain by 
the Bethe ansatz technique and to discuss the low temperature properties contributed by impurities 
with the use of the Landau-Luttinger description. The quantum inverse scattering method ]3lJ is 
adapted to solve the eigenvalue problem of the transfer matrix which determine the Hamiltonian 
of the impurity model. As is well- known, the scattering matrix for two electrons in the Hubbard 
chain is 

sin ki — sin kj — i Y 

where is the permutation operator of the electrons. The wave function of the chain in region 
< xqi < XQ2 < • • ■ < xqn < G — 1 has the form: 

^a u a 2 ,---^ N (xi,X 2 , • ■ • , X N ) 



N 

£ps r Aa Q1 ,..., aQN (rpQikpQi r --,rpQNkpQ N )exp[i'^2rp j kp j x j ]] (2) 

n---rjv— ±1 

where £_p = 1 (—1) when the parity of P is even (odd), e r — YljLi r in which r takes the value +lor 
— 1, which gives that A... aulJj ... (• • • fcj, fcj ■ • •) = Sij(kikj) x4... ffjif ,.... (• • • kj, fcj • • •) . The impurities 
are coupled to the two ends of the system and they are introduced by the scattering matrix with 
electrons, 

M,,,-**-*"-^', (3) 
sinfcj - iC L . R - % -f 

so that the boundary matrices take their form as 

n s sin kj — iCl — ^-Pl i 
R L (kj ,aj) = e^ L ^ > 3 - | r -^ , 4 

Rr (~kj, aj) = e -«MC+i) e Wfc,) sin k J ~ iCR ~ J Pr i ( 5 ) 
V 3 J> smkj+iC R +fP Rj 1 ; 

where the interacting parameters Cl,r is arbitrary. The boundary R matrix satisfies the reflecting 
Yang-Baxter equation and the eigenvalue problem of the system can be reduced to the form 

Tr T [T(\)T-i(-\)]\ x=s .. d> 



i - sin kj sin kj + iC'l + sin kj + iC R + ^ 
i^j — sin kj sin kj — iCl — ^jf sin kj — iC R — ^ 

.g-*{¥>i,(fei)+¥>H(fej)} e 2ifc i (G+l)^ 



with the transfer matrix defined by 

T (A) = S T j (A) 5V (A) • • • S T j-i (A) SV7+1 (A) • • • 5 r jv+i (A) , 

and S T i (A) = (A — sin A:; — iUP T i/2) / (A — sinfc; — iU/2) with sinfc = iCl, sinfcjy+i = iC R 
for I = 0, 1, 2, • • • , N + 1. Notice that the state function $ is related to the coefficient 
Aa Q1 ,a Q2 ,---,a QN (rpQ 1 kpQ 1 ,rpQ 2 kpQ2, ■ ■ ■ ,r PQN kpQ N ), which is dependent on also the spins of 
magnetic impurities and is suppressed for brevity. Pjo and PjN+i are the permutation operators 
between the magnetic impurities and the conduction electrons. The Hamiltonian of the model has 
the form 

H=- Y G to G 3° + U Y1 n A n H + ■ d L + Vim + JrSg ■ d R + V R n G (7) 

where dh and d R arc the spin operators of the impurities with the spins 1/2. The inter- 
acting constant Jl r and the scattering potential Vl r of the impurities can be expressed 
as J L . R = W L m/3 = -2/(1 + C 2 R ) when U =' 2 and J L . R = WV L , R /(A + U) = 

U (u 2 - AC 2 LR - 8 + ^Jl&Cl R - SC 2 R U 2 + MC 2 R + C/ 4 j / (8 - 2U 2 ) when U ^ 2 where we 

have assumed that the electrons scatter with the impurities with the small momcntums. From the 
Bethe ansatz equations of this system with the impurities, an immediate consequence is that the 
total momentum of the system has the finite size correction Y^j=i Sz=o l,r ® (9.7> ^ u + I $G) 
, contributed by the magnetic impurities coupled to the open boundary fermions system. For the 
Hubbard ring, Ij and J a are consecutive integers ( or half-odd integers ) centered around the origin 



and satisfying J2j kj = of the ground state. Now the magnetic impurities affect the distributions 
of Ij and J a and the index a runs from — M to M which is different from the periodic case. 

In the thermodynamical limits: N — > +00, G —* +00, M — > +00 with the ratios N/G, M/G 
kept finite, the real numbers k and A are distributed continuously in the ranges — Q to Q < tt 
and — B to B < +00 with the density functions p (k) and a (A) , respectively. The distributional 
functions satisfy the coupled integral equations 

9 ( k ) + Y) a (q, 2u + Q) = + ^ / a (A) [a {q-X,u) + a(q+ A, u)] dX, (8) 



ff ( A ) + ^E C <T{\')a{X-l\',2u)dX' = ^^- + \ J2 P P(k)a(\-lk,u)dk (9) 

with <j g (A) = a(X,u)+a (A, u + C{) + a (A, u — C/) + a (A, u + G#) + a (A, it — G#) where a (g, 77) = 
7r ~ 1? 7/ (<? 2 + ■ Q an d -B are the cutoffs of k and A , respectively, for the ground state and 
they satisfy that f_gp(k)dk = N/G and J B g a (A) dX = M/G . The ground state energy is 

E/G = -2 J Q Q p(k) cos kdk. For the half-filled band Q = it and AT/G = 1. The number of the 
down spins is a half of that for the conduction electrons in the system. To see this, by Fourior 
transformation of the second Bethe ansatz equation and assuming that the distributed functions 
are the even functions, we have that J_q p (k) dk = 2 a (A) dX with the choice of lu — and 
G — > +00. Then the magnetization is that S z = (N + 2 — 2M) /2 = 1, which is different from the 
pure system. The density functions are 

^ j 2cosh(ua;) 2G1+H(w,2«)' 1 ' 



1 cosfc f 00 , cos (wo) J (v) p G (k) 1 cos A: f°° , ct g (cj) e -«?"H«"l 

p(/c) = 1 / v ; y - - H / r — (11) 

PV ; 2tt tt J l + e 2u " 2G 2G 2ir J_ QC l + a(uj,2u) V ' 

where p (k) = cos k r a (?> + Ci) an d we have used that u > and Jo is the Bessel 

function. The terms with factor 1/(2G) describe the finite size corrections due to impurities. 
Then the finite size correction for the ground state energy due to the impurities is 

e' = r dk cos k P G (k) - r a / iM ^ h (12) 

J-* J -00 2wcosh(ww) 

where a G (ui) is the Fourior transformation of the function <j G (A). In order to investigate the 
conductivity properties of the ground state, we should computer the chemical potentials p+ and 
jit_ (see Ref. [22]). By the duality of the particles and holes in the system, for the case of the 
half-filled band, it should be derived that A = U — 2/i_ with 

^r = -*f dk cos k[p{k)-p (*;)'] + I (13) 

where the density function p(k)' is determined by the coupled integral equations (||) and 

dX'a (A') a (A - A', 2u) + a (A) = / dkp (k) a (A - q, u) + — ^ ~ 2& ^' ^ 



2G 

It is due to that a hole with kh — n and A — > 00 appeared comparing with the 'half-filled' band and 
this gives the term a (A, u) /G in the above equation and term 2/G in relation (113). The interesting 



A 



thing is that the impurities do not change the conductivity property of the system for the ground 
state although the density functions have been corrected into the forms (JlOj) and ([n]). 

The other physical properties of the system can be discussed by evaluating its thermodynamics 
from the Bethe ansatz equations. As an alternative, both more practically and more physically, 
we now use the picture of the Landau-Luttinger liquid, putting forward by Carmelo and co- work- 
ers |52| , to determine the magnetic contributions to the specific heat, susceptibility and Kondo 
temperature. The densities of the states of the quasiparticles at the Fermi energy Ep = are 
described by 



iV c (0) = 



1 



2 r KV r 



1 



1 Sp(ko) 



2G p(k ) 



N s (0) 



1 



2irv s 



1 



1 5a (00) 
2G a (00) 



where v c and v s are the velocities of the charge and spin fluctuations, respectively, p and a are 
the distributed functions in the infinite limit. And dp, 8a are the finite-size corrections due to 
impurities. It means that the density functions are denoted by p + 5p/ (2G) and a + 5a/ (2G) 
of parameters k and A , respectively. The densities of the states determine the low- temperature 
behaviors of the system by the use of the standard expressions of the Fermi-liquid theory and we 
get the contributions to the Kondo effects due to the magnetic impurities as the following form, 



5C 



irSp(k ) 



T 



ir5a (00) 



6Gv c p (fco) 6Gv s a (00) 



T. 



(14) 



where the specific heat comes from the charges and the spins. The finite-size correction to the 
susceptibility is 8\ = XoSa (°°) / [Ga (00)] where xo is the susceptibility in the bulk. The Kondo 
temperature Xfc, corresponding to the Fermi temperature in the local Fermi liquid generated by 
the impurities, can be derived from the impurity specific heat contributed by the spin sector. It 
has the form T k = 2Ga (00) / [irnSa (00)] where n is the density of the electrons in the system. 
Furthermore, under the limits of the parameters Cr = Cl = 0, the specific heat due to the spin 
of the impurity, the susceptibility and the Kondo temperature are 



sc 



5ttT 



6Gv s Jq («7r) 



.Sx = 



GJ (*7T) 



Xo,T k 



2GJ (wt) 
5mr 



In this case we have that Jl,r — 4Z7Vl,h/(4 + U) = —U. Notice that the Kondo temperature Tk is 
linear in the 1 jn and this property is similar as the one in Ref. [29] . The specific heat contributed 
by the charge of the impurities have the form: 



6C r .= 



7T [In 2 + 2(3 (3/2)] 



6Gv r 



2 J2Zi (! + 4« 2 Z 2 r 1/2 - 1 



T 



(15) 



with Cr = Cl = U/4 and under the other cases the specific heat SC C due to impurities have 
the similar expressions at the low temperature. Here the (3 function is defined by f3 (x) = 

EfcLo(-l)* /(* + *)> and (3/2) « 0.429 2. 

In summary, we have studied the low energy properties of the Kondo problem in a ID chain 
of strongly interacting electrons described by the open boundary Hubbard model. The magnetic 
impurities are coupled to the ends of the system and they are introduced by the scattering matrices 
with the electrons. This is the first approach to study the impurity properties in the open boundary 
Hubbard chain. The finite size correction of the ground state energy is obtained and we find 
that the magnetic impurities do not change the conductivity for the ground state although the 
distributed functions have been changed. By the use of the Landau-Luttinger liquid description 
for the Kondo problem, we obtained the expressions of the low-temperature specific heat, the 
susceptibility and the Kondo temperature. The Kondo temperature is inversely proportional to the 
density of electrons. Finally, we point out that the integrability of the impurity model is preserved 
due to the boundary matrices satisfy the reflection equation. The corresponding Hamiltonian 
has been written down explicitly when the electrons scatter with the impurities with the small 



momcntums, which has a simple and compact form. The general expression of the Hamiltonian 
can be constructed with the use of the quantum inverse scattering method and we wish to remain 
it as an open problem for the further investigations. 
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